CDM and Baryons as Distinct Fluids in a Linear Approximation for the Growth of 
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A single fluid approximation which treats perturbations in baryons and dark matter as equal 
has sometimes been used to calculate the growth of linear matter density perturbations in the 
Universe. We demonstrate that properly accounting for the separate growth of baryon and dark 
matter fluctuations can change some predictions of structure formation in the linear domain in a way 
that can alter conclusions about the consistency between predictions and observations for ACDM 
models versus modified gravity scenarios . Our results may also be useful for 21cm tomography 
constraints on alternative cosmological models for the formation of large scale structure. 
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I. INTRODUCTION 

Since the discovery that the universe is apparently 
dark energy-dominated causing an observed acceleration 
([H, [|J, also inferred indirectly on the basis of other ob- 
servational constraints, i.e. see, for example Q), a vast 
expenditure of effort has been made towards possible ex- 
planations of the acceleration. The standard paradigm of 
cold dark matter with a cosmological constant (ACDM) 
in the cosmological framework of general relativity (GR) 
accommodates all experimental evidence, and remains 
the simplest and most economical cosmological model 
consistent with the data. Though the ACDM cosmology 
fits all the present data, issues such as the hierarchy and 
coincidence problems remain which highlight the issue of 
how the acceleration can be realized in a fully consistent 
theoretical framework. This has led to a consideration 
of alternative cosmological models, including models in 
which gravity varies away from GR on large scales. As 
observations become increasingly precise, the ACDM pic- 
ture will be put through ever more rigorous tests in the 
effort to constrain new physics. It is important there- 
fore to have accurate theoretical frameworks by which to 
judge whether observations may indicate a discrepancy 
with the predictions of the standard model. 

By now a standard way to constrain various cosmo- 
logical alternatives is via an exploration of the growth of 
linear matter perturbations for various redshifts. These 
perturbations have been parameterized via a growth in- 
dex (see for example ) , which has a specific value for 
ACDM (to first order in deviations from a purely CDM 
dominated universe, where the deviation is due to a cos- 
mological constant, this index has been estimated to be 
simply 6/11). 

The growth is typically found in the following manner. 
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One defines a matter overdensity given in fc-space by 6 = 
5p/ p with p being the background matter density. Using 
the standard Einstein equations, one finds the dynamical 
equation typically called the growth equation, for a single 
component matter field in a matter dominated universe 



5 + 2H5-A-kG nP 5 = Q 



(1) 



where the overdot is a derivative with respect to coordi- 
nate time and we have dropped the k index. This relation 
is given in the synchronous gauge where it holds on all 
scales (for gauge related issues see for example [Tol - [l7l |). 
From here one can define a function g = dlnJ/dlna = 
^m(a) 7 which leads to the equation 



H 



(2) 



where the prime denotes a derivate with respect to the 
natural log of the scale factor. The function g can be 
identified as the growth factor for matter density pertur- 
bations. The solution to this for a flat universe with a 
dark energy equation of state w\ is given, to first order 
in the expansion parameter 1 — Q,\ by 
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3(^a - 1) 
6w\ — 5 



(3) 



This relation reduces to 7 = 6/11 for the case of ACDM. 
The growth factor g defined in (2) is affected only by 
CDM overdensities. In our paper we will consider the 
effect of considering both CDM and baryonic perturba- 
tions on the growth factor. 

Therefore we focus on two facts: 

• The matter content of the universe is not solely 
composed of cold dark matter, as the baryonic con- 
tent is roughly one- fifth that of dark matter [l8j . 
If one writes separate growth equations for dark 
matter and baryonic matter densities, each contain 
a source term involving the gravitational potential 
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which is a function of the full matter content, in- 
cluding both baryons and dark matter (for exam- 
ple, the third term in Eq. (JlJ arises from the Poisson 
equation for the gravitational potential). 

• Because baryonic matter has a non-zero sound 
speed, dark matter and baryonic matter perturba- 
tions obey different dynamical equations |19l — 1 211 ] - 

Here we explore the consequences of properly incorpo- 
rating both of these effects, and quantify and compare 
differences in the perturbativc densities when the full set 
of baryonic plus dark matter equations are solved versus 
the case when baryons are ignored. One of the central 
purposes of this effort is to compare the relative differ- 
ences between results of these two approaches compared 
to the differences obtained when different cosmological 
models are explored, in order to determine the sensitiv- 
ity to cosmological model dependence versus the need to 
properly account for baryons. 

To explore these effects we perform calculations under 
both the standard ACDM scenario and a modified grav- 
ity model, for which we choose the DGP [22| model. We 
compare bias factors, total matter density perturbations 
and the growth factor in these to cosmological scenarios 
scenarios to explore the sensitivity to not including bary- 
onic perturbations. We find that an accurate treatment 
of baryonic fluctuations will alter quantities like the bias 
factor and the total matter density fluctuation even in 
the linear regime in a way that can exceed the change 
induced in the quantities by varying the the background 
cosmology. On the other hand we find that the growth 
factor, defined as §^ is relatively insensitive to inclusion 
or non-inclusion of baryonic dynamics over linear scales. 
It is on the other hand very sensitive to the background 
cosmological model. 

In Section II we briefly present the formalism for the 
full set of dynamical equations. In Section III we present 
our results, and finally in Sec. IV we conclude with a 
brief discussion of their implications. 



II. CALCULATION 

The full set of coupled linear differential equations 
for the growth of perturbations in dark matter (S c ) and 
baryons (S b ) along with radiation (5 ra d) is 



S c + 2H(z)S c 
S b + 2H{z)5\ + c 2 s k 2 

&rad 



H 2 (f c S c + f b S b + f r S rad ) 

H 2 {f c S c + fbfa + frSrad) 



(4) 



f c = fl c (z) fb = Ctb(z) 

fr = tlrad(z),0 (z > Z dec , Z < Z dec ) (5) 



where tt c , tt b and VL rad arc CDM, baryons and radiation 

energy density at a given epoch. We define c s = (^j 

[20j . where the subscript S stands indicates that c s is 
defined at constant entropy, S. At low redshifts c s is pri- 
marily due to baryonic pressure and at high-rcdshift (prc- 
decoupling redshifts), the contribution is mainly from ra- 
diation pressure. At high redshifts we use Eqns.(|3]) and 
([5]) to solve for the baryonic and dark matter perturba- 
tions. We evolve the perturbations in each Fourier mode 
fc, starting from the epoch of horizon entry for that par- 
ticular mode. Therefore our boundary values for S b and 
5 C , corresponding to a mode k, are set at the epoch of 
horizon entry, z en t , such that 
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fco 



SScOBE 

a{z eq )H(z ent ) 



Sc,b(z en t,k) = A 



A 



(fc < k e q) 



(k > k e 



5 c ,b( Z ent,k) 



S c ,b(z ent , k)H(z ent ) (z ent > z eq ) (6) 

S c ,b(Z en t, k)a(z en t)H (Zent) (z ent < z eq ) 



In Eq.©, modes that enters the horizon at the epoch of 
matter and radiation equality (z eq ) are denoted by fco- 
Fluctuations are normalized by using the temperature 
fluctuation over angular scales of 7 degrees at the surface 
of last scattering measured by the COBE mission which 
is denoted by 6cobe |23[. The spectral index of the 



primordial fluctuations (coming from very early times) is 
set by n s over all Fourier modes. 

The boundary conditions for S c ^ b described by Eq.© 
come from the following argument. Given a primor- 
dial power spectra of shape P,(fc) cx fc" s , fluctuations 
grow during the radiation dominated epoch such that 
5 oc a 2 . Therefore when a given mode enters the horizon, 
its power is described by P ent (fc) oc a 2 Pi(k) ~ fc" s ~ 4 if 
Zent > z eq [U, and we replace 5 c ,b{z ent , fc) oc \JP en t{k). 
Note that if n s « 1, then fc 3 P eTlt (fc) is constant at hori- 
zon crossing. Similarly, we find that, P en t{k) oc fc" 3 cor- 
responding to the modes which enter the horizon after 
matter-radiation equality. 

The calculation of c s involves matter and radiation 
temperatures along with their fluctuations. We calcu- 
late the matter temperature, T ma t(z), at a given epoch z 
using the following equation [25j . 



dT„ 



df 



+ 



2H(z)f mat 
x e (t) 



(7) 
(8) 



The radiation temperature, T(z), at a given epoch z is 
estimated by the standard relation T(z) = Tq(1 + z), with 
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To given by the present CMB temperature [26[ . Fluctua- 
tions in the matter temperature after mechanical decou- 
pling (z > 1100), St, are calculated using 



cISt 



2dS b 
3~dt 



x e (t) 



-^s 7 

T 



(9) 



fDGP model we use il c h 2 = 0.12,n b h 2 = 0.023, h = 0.69. 
Eqn.(4-5) of [29[ are incorporated to compute the mod- 
ified expansion rate Hdgp{ z )- We have plotted up to 
z ~ 1 to highlight the effect at low redshift, where one 
would expect modifications of gravity designed to mimic 
dark energy to be most relevant. 



St thus computed is in turn used for the calculation of 
the sound speed post recombination, and therefore the 
modified baryonic growth equation becomes [27j | 



2H(z)S b 



\H 2 {f c S c + f b S b ) 
-(S b + S T ) 



a 2 /i 



(10) 



For scales which enter the horizon before the epoch of 
recombination, z rec , matter temperature fluctuations are 
described by St = St^ at z = z ent . 

As we would like to examine growth not only in the 
standard cosmology, but in a modified gravity scenario as 
well, we will now discuss how the calculation needs to be 
altered. In the DGP scenario, the CDM and baryon per- 
turbation equations are modified such that in the source 
term on the right hand side of Eq.Q, the factor H(z) 2 
is replaced by Hp GP g(a, k). In this case Hdqp(z) is the 
modified background expansion rate and g(a, k) is the 
factor by which Newton's constant gets modified under 
the new gravity scenario. 

We use [28[ and [29| to construct Hdgp( z ) and g(a, k). 
Using these modifications and Eqns.(|3][S]), we calculate 
the growth of perturbations in the DGP theory up to a 
scale corresponding to k < 0.05Mpc~ x . We chose to re- 
strict ourselves to these scales in order to avoid complica- 
tions due to non-linear PPF parameters as described in 
J28|. The matter temperature perturbations in Eqns.(0- 
riTJj) also get modified by replacing H(z) with Hdgp(z). 

The input parameters of our calculation are n s , St~ , 
fi^c h, fife, r^tot and z eq . Additionally, the epoch 
of decoupling, Zd ec , is determined such that the pho- 
ton mean-free path is larger than Hubble distance, or 
A 7 = - 1 ~ cH^ 1 . For simplicity, we set tt u = 0, 

= and allow a sharp drop in optical depth of pho- 
tons at the epoch of mechanical decoupling. In the next 
section we will describe our results using WMAP9 values 
[26|, along with additional DGP fits. 



III. RESULTS 

In Figure [T] we represent the background expansion 
rate with respect to redshift, H(z), in different cosmologi- 
cal cases. We consider fCDM, fDGP (flat ACDM and flat 
DGP), and oDGP (open DGP) models. For cosmological 
parameters we use Q c h 2 = 0.12, fl b h 2 = 0.023, h = 0.69 
from WMAP9, and apply those to both fCDM and fDGP 
models. For the oDGP models we use fl c h 2 = 0.099, 
flhti 2 = 0.023, fl k = 0.03 and h = 0.76 [H, and for the 




FIG. 1: Background expansion rates (H/Ho) as a function of 
redshift based on different fCDM (flat ACDM), fDGP (flat 
DGP) and oDGP(open DGP) cosmologies. Relevant cosmo- 
logical parameters are described in the text. 

In Figure[2]we present the bias, defined as b(z, k) = j 1 -, 
with respect to redshift. We choose three length scales 
corresponding to Fourier modes k — 0.005, k = 0.01 and 
k = 0.05 in units of Mpc -1 , where b(z, k) is represented 
in those regimes respectively by solid, dotted and dashed 
lines. We choose fCDM, fDGP and oDGP cosmologies 
described by the same cosmological parameters as in Fig- 
ure [T] For a single fluid model the bias is unity by con- 
struction. From Figure[2]it is evident that in all scales ex- 
plored, the difference in bias between the DGP and CDM 
models is less than their deviation from unity. This de- 
viation (with respect to unit bias in a single-fluid model) 
increases both with redshift and diminishing scale. For 
low redshift (z - 0.2) and k = O.OlMpc -1 , the differ- 
ence between the bias calculated from fCDM and fDGP 
models (using WMAP9 parameters) is about 0.04% with 
a difference of 4.2% in the background expansion rate. 
The bias becomes close to 1% near z ~ 1, while the dif- 
ference in bias between fCDM and fDGP (using WMAP9 
parameters) is only up to 0.02% for k = 0.05Mpc _1 . 

The bias can be directly related to the total den- 
sity fluctuation S = f c S c + f b S b such that S = 
( fc + b(z, k)f b ) 5 C , where b{z, k) represents the bias at a 
given epoch and a given scale k. In Figure [3] we display 
Stot{z,k) as a function of redshift in the linear regime, 
k = 0.005Mpc _1 . The same normalization was used for 
all cosmologies, set by Eq.©. We use dotted lines to re- 
fer to the single fluid (s-f) models and solid lines for the 
baryon+CDM fluid models. We note that at low red- 
shift, the difference in the oDGP and fCDM models is 
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FIG. 2: The bias as a function of redshift in different cos- 
mologies for scales corresponding to k — 0.005(solid lines), 
k = 0.01(dotted lines) and k = 0.05(dashed lines) in units of 
Mpc" 1 . 



comparable to that between the fCDM, s-f and two fluid 
models. This is an intriguing conclusion which suggests 
the importance of the two-fluid treatment in order to cor- 
rectly use structure formation observations to constrain 
cosmological models. The significance of 8 tot is that it is 
a scale and cosmology dependent quantity. Observation- 
ally, future weak lensing surveys can estimate Stot but 
the interpretation of observational data must be made 
by properly incorporating bias on all scales of interest. 
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FIG. 3: The total matter fluctuation S = f c S c + fbSb is shown 
as a function of redshift under different cosmologies at a scale 
corresponding to k = 0.005Mpc _1 . Solid lines indicate two- 
fluid models and dotted lines represent single fluid (s-f) mod- 
els. 



In Figure 21 we plot growth factor only due to CDM 
c (solid lines) with respect to redshift for a cho- 

1 . Note that 



" c dlna 

sen scale corresponding to k 



0.005Mpc 

we over plot (dashed line) the parameterization D,m f30| 
which agrees quite well with the numerical estimates for 
fCDM. We therefore conclude that g c is roughly scale- 



independent, but is sensitive to the background cosmol- 
ogy 




FIG. 4: The growth factor due to only CDM described by 
g c — is plotted as a function of redshift for differ- 

ent fCDM°(flat ACDM), fDGP (flat DGP) and oDGP(open 
DGP) cosmologies. This is done at a single scale k = 
0.005Mpc~ .Two component models are presented in solid 
lines and the dash-dotted line referred as analytic indicates 
f2 m (z) 7 such that 7 = £ for fCDM using WMAP9 parame- 
ters. 



We can also consider the evolution of the total growth 



factor due to CDM and baryons as g t0 



dlnSto 
dlna 



with re- 



spect to redshift at a scale of k = 0.005Mpc 1 , to check 
to see if there is any difference in using this value instead 
of the CDM growth factor. Wc define 8 to t = S °f'+^ fb ■ 
We find that both parameters g c and gtot are weakly de- 
pendent on the selection of scale k., with g to t to be com- 
paratively more rigid over a range of various k values. 
However in general g c ~ gtot so there is no significant 
handicap in using g c to constain cosmological models. 



IV. DISCUSSION AND FUTURE DIRECTIONS 

In this work we have examined the impact on calcula- 
tions of the growth of structure through the use of a sim- 
ple single matter fluid approximation vs. a model that 
correctly incorporates baryons and cold dark matter in a 
two component analysis. We have then compared the dif- 
ference in matter perturbations in the standard ACDM 
cosmology and modified DGP gravity. 

Among the quantities wc discussed in our paper, the 
growth factor g = 4jr^ is measured observationally from 
galaxy redshift surveys [3l|. We found the growth factor, 
g to be rigid with variation in k under both ACDM and 
DGP cosmologies, and largely independent of baryonic 
dynamics. 

Using weak lensing to measure the total 8 = f c 8 c + fb8b 
is a way to get a handle on total matter density fluctua- 
tions but this quantity is sensitive to baryonic dynamics 
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Q and depends on the relevant scale of structure forma- 
tion even in the linear regime. In addition, we find that 
at low-redshift (0.7 < z < 1.5) the effect of baryonic dy- 
namics can be comparable to that introduced by modify- 
ing the underlying cosmology. We show in Figure [3] that 
for z ~ 1, modifications due to two component model 
is comparable to that due to modifying the cosmology. 
Therefore accurate inclusion of the baryonic dynamics 
is required for interpreting observations associated with 
this quantity 

It is also worth noting various studies of the growth of 
large scale structures in different modified gravity con- 
texts have been performed using a single growth equa- 
tion (for a sample see [32T - I4]] ] ) . Preliminary results, to 
be described in a future work suggest that the use of a 
proper two component fluid formalism can significantly 
weaken the ability to distinguish between cosmological 
predictions in such models. 

Finally, another area of cosmology which has attracted 



a great deal of interest recently, and relies on accurate 
calculations of the evolution of density perturbations, is 
the signal arising from the 21cm spin-flip transition of 
neutral hydrogen (for recent reviews, see j42T - |i3 |). One 
may study, for example, the perturbations of the bright- 
ness temperature of the CMB over a large redshift range 
in the so-called Dark Ages, which rely on the density per- 
turbations of hydrogen. These perturbations are seeded 
by dark matter, and therefore a precise calculation merits 
the inclusion of the full baryonic plus dark matter system. 
We are currently using the formalism we have described 
here to investigate how it will impact upon conclusions 
one may draw from the use of such observations. 
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